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ABSTRACT 

We present high signal precision optical reflectance spectra of the large Kuiper 
belt objects 2005 FY9 and 2003 EL61. The spectrum of 2005 FY9 exhibits strong 
CH 4 -ice bands. A comparison between the spectrum and a Hapke model indicates 
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the CH4 bands are shifted 3.25 ± 2.25 A relative to pure CH4-ice, suggesting the 
presence of another ice component on the surface of 2005 FY9, possibly N 2 -ice, 
CO-ice, or Ar. The spectrum of 2003 EL61 is remarkably featureless. There is a 
hint of an 02-ice band at 5773 A; however, this feature needs to be confirmed by 
future spectroscopic observations of 2003 EL61 with a higher continuum signal 
precision, sufficient to detect a second weaker 02-ice band at 6275 A. 

Subject headings: Kuiper Belt — techniques: spectroscopic 



Introduction 



Measuring the surface composition of Kuiper belt objects (KBOs) may provide clues 
about the composition and environment of the primordial solar nebula as well as the im- 
portant evolutionary processes occurring in the outer Solar System over the last 4.5 Gyr. 
Unfortunately, there are only a handful of KBOs that are known to ex hibit ice absorptio n 
bands in their spectra. H^O-ic e bands are seen in spectra of 19 96 TO66 (IBrown et al.lll9 99). 



Varuna (jLicandro et al.l 120011 ) , Quaoar (IJewitt fc Luu 



1980 



2004 a|). CH_4-ice bands are seen in spectra of Pluto (ICruikshank et al.l 1 19761 ; iFink et al. 



Grundy fc Fink 



2004), and Orcus (Fornasier et al. 



19961 ). N eptune's satellite T riton, which may be a captured KB O 



Cruikshank et all Il993h . Eris (brown et all l2005f ). and 2005 FY9 JLicandro et all l2006h 
. There are perhaps a dozen objects that exhibit spectra with no ice absorption bands 
( IDoressoundiram et al.l 120031 ; iFornasier et al.l l2004bl ) 



The recent discovery of extraordinarily bright and large KBOs opens up a new opportu- 
nity for the physical and chemical studies of KBO surfaces with high signal precision optical 
reflectance spectroscopy. Bright objects make it possible to either carry out in-depth physi- 
cal modeling of ice and mineralogical absorption bands or set stringent upper limits on the 
presence of such bands. 

Here we present high signal precision optical reflectance spectra of KBOs 2005 FY9 
and 2003 EL61. Both objects are candidates for membership in the newly defined class 
"dwarf planet." The diameter of 2005 FY9 is estimated at 1600 km (IBrown et al.l 120061 ) 
while 2003 EL61 appea rs to be a highly elongated ellipsoid with axes of 1950 x 2500 km 
( IRabinowitz et al.ll2006l ). For comparison, Pluto has a diameter of 2350 km. 2005 FY9 
and 2003 EL61 are among "scattered-near" objects (i.e. non-resonant, non-planet-crossing 
objects with Tisserand para meters less than 3 , relative to Neptune) in the Deep Ecliptic 
Survey classification system (jElliot et al.ll2005l ). 2005 FY9 has a perihelion distance, q, of 
38.6 AU, an aphelion distance, Q, of 52 AU, a semi-major axis, a, of 45.3 AU, an inclination 
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angle to the ecliptic, i, of 29°, and an eccentricity of 0.15. 2003 EL61 has q = 35.1 AU, Q 
= 51 AU, a = 43.1 AU, i = 28°, and e = 0.18 For comparison, Pluto has q = 29.7 AU, Q = 
49 AU, a = 39.5 AU, i = 17°, and e = 0.25. At the time of observation, 2005 FY9 and 2003 
EL61 were both near aphelion. 



2. Observations 

We obtained spectra of 2005 FY9 and 2003 EL61 on 2006 March 5 UT with the 6.5- 
meter MMT telescope on Mt. Hopkins, AZ, the Red Channel Spectrograph, and a 1200 x 
800 CCD. A 150 g mm -1 grating and a 1 arc sec slit width provided wavelength coverage of 
5000 to 9500 A in first order, a dispersion of 6.38 A pixel -1 , and a fwhm resolution of 20.0 
A. A LP495 blocking filter eliminated contamination from the second order. 

There were high, thin cirrus clouds through most of the night and the seeing was ~ 0.8 
arc sec. The KBOs were placed at the center of the slit and tracked at KBO rates. 

We used the Im age Reduction and Analysis Facility (IRAF) and standard procedures 
( IMassey et al. Il992h to calibrate and extract one-dimensional spectra from the two-dimensional 



spectral images. Specifically, the electronic bias of each image was removed by subtracting 
its overscan as well as a bias picture. Pixel to pixel sensitivity variations were removed 
from each image by dividing by a normalized twilight flatfield image. Extraction of one- 
dimensional spectra from the images was done with the apall task in IRAF. HeNeAr spectra 
were used to correct for flexure and obtain an accurate wavelength calibration. Our wave- 
lengths are accurate to ~ one-tenth of a pixel or ~ 0.7 A. We removed telluric b ands and 



Fraun hofer lines from the KBO spectra by observing the solar analog HD 112257 (jHardorp 



19821 ) at airmasses very close to 2005 FY9 and 2003 EL61, and then then dividing the KBO 
spectra by the normalized solar analog spectra. Typically, the airmass difference between 
the KBOs and the solar analog was < 0.05. The 2005 FY9 and 2003 EL61 spectra in the 
analysis below are an average of four 10-minute and five 10-minute exposures, respectively. 



Results 



3.1. 2005 FY9 



We plot our reflectance spectrum of 2005 FY9 in Figure 1 (bottom). The tick marks 
correspond to previously reported absorpt ion maxima of CILt-ic e (Table 1). The middle and 
top spectra in Figure 1 are of 2005 FY9 (ILicandro et al.ll2006l ) and Pluto (I Grundy fc Fink 



-4- 



1996). All spectra a re normalized to 1 at 6500 A. For the purpose of comparison, the 
Licandro et al.l (120061 ) 2005 FY9 and Pluto spectra are offset by 0.4 and 0.8, respectively. 
Clearly, 2005 FY9 exhibits much deeper CH 4 -ice bands than Pluto. 



Our spectrum differs from the iLicandro et al.l (120061 ) spectrum in two ways. First, our 
spectrum exhibits additional subtle absorption bands at 5400, 5800, and 6000 A. These 
bands have not been seen before in astronomical or laboratory spectra of CH 4 -ice. Because 
the bands are close to the wavelengths of gas phase CH 4 bands and because there are strong 
CH 4 -ice bands at longer wavelengths in the spectrum of 2005 FY9, we think the new bands 
are due to CH 4 -ice . In a ddition, it appears the band at 6200 A is CH 4 -ice as well. It is in 
the ILicandro et al.l (120061 ) spectrum, but they do not identify it a s CH 4 -ice. Second, we find 
a continuum slope between 5500 and 6500 A of 10%/1000A while ILicandro et al.l (120061 ) find 
a slope of 13%/1000A, with both spectra nor malized to 1 at 6 500 A . We used HD 112257, 
a G2V star, for our solar analog star whereas ILicandro et al.l (120061 ) used BS 4486, a GOV 
star, for their solar analog star. An examination of the V— R colors for GOV and G2V stars 
suggests the slope difference is due to the use of different solar analog stars. 

In order to place constraints on the CH 4 grain size as well as the presence of any 
additional ice components on the surface of 2005 FY9, we calculated m odel CH 4 -ice spectra . 
Our models use laboratory optical constants for pure CH 4 -ice at 30 K (IGrundy et al.ll2002l ) 
and arbitrary optical constants that absorb more at blue than red wavelengths, and thereby 
reproduce the observed reddish slope for the 2005 FY9 continuum. The ice and reddening 
agent were mixed in the model at the molecular level (i.e. by means of a weighted ave r age o f 
the ice and reddening agent optical constants). We used Hapke theory (jHapkd Il98ll . Il993l ) 
to transform the optical constants into reflectance spectra. We used Hapke parameters h 
=0.1, B G = 0.8, 9 = 30°, and P(g) = a two component Henyey-Greenstein function with 
80% in the forward scattering lobe and 20 % in the back scattering lobe and both lobes 
having asymmetry parameter 0.63. These values are comparable to numbers used previously 
to model the surface of Pluto (IGrundy fc Buidl200ll ). 



Figure 2 compares the spectrum of 2005 FY9 (black line) and a Hapke model with a 
grain size of 0.55 cm (red line). Although the model fits the core of the 8897 A band, it has 
too little absorption at all the other bands. Next, we fit the core of the weaker 7296 A band 
with a 1.8 cm grain size. Then, we fit the core of the even weaker 8000 A complex with a 
2.5 cm grain size. These models illustrate a trend, weaker bands (which see deeper into the 
surface) aways seem to need bigger grain sizes, or equivalently, greater optical path lengths 
in CH 4 -ice. 



From the above discussion, it is clear a single grain size will not fit the observations. In 
Figure 3, we present a Hapke model that uses two grain sizes, 6 cm (97 % by volume) and 
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0.1 cm (0.03 % by volume), in an intimate mixture. The inability of the two-grain size model 
to fit the CH 4 -ice band at wavelengths less than 7000 A in Figure 3 is due to the lack of 
laboratory optical constants. We note that our two-grain size model is not a unique solution. 
Although this model fits the data, other models (such as courser grains underneath a coating 
of finer grains) would probably fit the data just as well. In addition, it is important to realize 
that grain sizes larger than 1 cm are probably a measure of the spacing of fractures or voids 
in the solid surface rather than indicating a surface covered by golf-ball sized particles. 

Although we can fit the shape and depth of the absorptions, there is a small but signif- 
icant difference between our models and the telescope data. Specifically, the maxima of the 
CH 4 absorption bands in the spectrum of 2005 FY9 reside blueward of the maxima in the 
pure CH 4 Hapke models. In Figure 4, we present a small portion of Figure 3 centered on the 
7296 A band to illustrate the shift. Such a shift is important because laboratory experiments 
show that the 8897, 8968, and 9019 A b ands in spectra of pure C H 4 -ice are shifted blueward 



by 27 A in a CH 4 /N2 = 0.8 % mixture (jQuirico fc Schmittlll997l ). Another experiment fi nds 



id 



the 8897 A band shifts ~ 17 A in a CH 4 /N 2 = 20 % ice mixture (IGrundv fc FinHll993 ). 



A cross correlation experiment provides a way to quantify the apparent shift in Figure 
3. Specifically, we wrote a Fortran program to shift the model spectrum from — 25 to +25 
A in 1 A steps. For each shift, the program finds the difference between the data and model, 
i.e. it calculates (y^i — y m ,i) where y^i and y m j represent the ordinate values of the data and 
the model spectra at wavelength i. Then, the program sums the squares of the difference 
over all N wavelength points between 7000 and 9300 A. In other words, we calculate 



N 

Ashift = ^2(yd,i - y m ,i) 2 - (i) 

i 

In Figure 5, we present a plot of the sums of the squares of the differences, i.e. A shi f t , as 
function of shift. We find a well defined minimum at 3.25 A. 

What is the uncertainty in the 3.25 A measurement? We note that HeNeAr spectra 
enable us to calibrate the wavelengths in the 2005 FY9 spectrum to an uncertainty of ~ 
1/10 of a pixel or ~ 0.7 A. Specifically, we find the average difference between the centroids 
of airglow lines in our data and the corresponding laboratory values is 0.7 A. A larger source 
of uncertainty comes from the noise in the spectra and the broadness of the cross correlation 
minimum in Figure 5. Propagation of the noise in the data and model spectra for the 3 A 
shift through equation (1) gives an ordinate error bar for the minimum point (3.25, 2.218) 
in Figure 5 of 0.042. As a result, the uncertainty in our shift measurement is defined by 
A shift < 2.218 + 0.042 = 2.260 (i.e. the section of the curve below the dashed line in Figure 
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5), corresponding to shifts between 1 and 5.5 A. In short, we find the Cr^-ice bands in the 
spectrum of 2005 FY9 are blueshifted relative to the model Hapke spectrum by 3.25 ± 2.25 

A. 



Licandro et al.l (120061 ) gave blueshifts of 2 - 6 A for CH4 bands in their Table 1, but they 



questioned the reality of the shifts because they are so small. Our independent measurement 
of shifts with sim ilar magnitude bolster s the case for the reality of the small shifts. In fact, 
when we take the iLicandro et al.l (120061 ) data and do a cross correlation experiment between 
our spectrum and their spectrum, we find no measurable shift between the spectra; but when 
we perform a cross correlation experiment between their spectrum and the Hapke model, we 
find their CLLt- i ce ba nds are blue shifted 2 A relative to the model. Our spectrum and the 
Licandro et al.l (120061 ) spectrum are consistent. 



3.2. 2003 EL61 



We plot our reflectance spectrum of 2003 EL61 in Figure 6. The continuum has a 
slope of essentially zero , consistent with B— V = 0.63 ± 0.03 and V— R = 0.34 ± 0.02 
(IRabinowitz et al.ll2006l ). Evidently, tholins are not present on the surface of 2003 EL61 
at an abundance seen on 2005 FY9. (Recall that both objects were observed on the same 
night and compared against the same solar analog spectrum.) We find no evidence of CH 4 - 
ice bands, but we can set an upper limit on the thickness of a global glaze of CH4-ice by 
applying Beer's Law, 



-2a x t 



(2) 



to the continuum near the 8897 A band. In the above equ ation, t is the thick ness of the 
glaze, aA are the Lambert absorption coefficients for CH4-ice (IGrundy et al.ll2002l ). the factor 
of two comes about because light passes through the glaze once on the way in and once on 
the way back out, (I/I )a is the ratio of light leaving the glaze to the light incident on the 
glaze. For (I/I )a=8897 = 0.97, we find t < 0.3 mm. 



H^O-ice bands are seen in near-infrared spectra of 2003 EL61 (IBarkume et al.l 120061 ) . 
Therefore, it is r easonable to look for evidence of solid-state photolytic or radiolytic chemistry 
(IJohnsonl Il990l ; IJohnson fc Quickendenl Il997l ; iDelitsky fc Land Il997l ) by seeing if 2 -ice is 
present on the surface of 2003 EL61. In Figure 6, we superim pose the 5773 and 6275 A 



02-ice bands in spectra of Ganymede (jCalvin fc Spencer 



1993) on top of our 2003 EL61 



spectrum. There is a tantalizing dip at 5773 A in the spectrum of 2003 EL61, but no dip at 
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the weaker 6275 A band. The continuum signal precision at the position of these bands is 
about 150 and so if the 5773 A band is real, it would take a continuum signal precision of 
about 300 to detect the 6275 A band. 

We estimate an upper limit to the thickness of a global glaze of 2 -ice by assuming the 
possible feature at 5773 A is real and then applying Beer's Law as we did for CH 4 -ice. For 
(I/Io)a=5573 = 0.99 and a.\=5j73 = 5.5 cm -1 , we find t < 0.01 mm. Our value for aA=5773 
comes fro m an integrated abs orption coefficient of 1.1 x 10 3 cm -2 and a band fwhm of ~ 
200 cm" 1 Landau et alill962h . 



4. Conclusions 

Ou r spectrum of 2 05 FY 9 exhibits strong CH 4 -ice bands in agreement with the spec- 



trum of iLicandro et al.l (120061 ). From a comparison of our spectrum and a Hapke model, we 
find the 2005 FY9 CH 4 -ice bands are shifted 3.25 ± 2.25 A blueward relative to the positions 
of pure CH 4 -ice bands. The shift could be due to the presence of another ice-component, 
possibly N2-ice, CO-ice, or Ar. Future higher resolution spectra of the individual CH 4 bands 
should determine if the weaker bands which penetrate deeper into the surface exhibit dif- 
ferent shifts than the stronger bands which do not penetrate as deep into the surface. Such 
observations could provide a technique to measure the CH 4 concentration relative to the 
other ice component as a function of depth below the surface of 2005 FY9. In addition, 
it is highly desirable to obtain a measure of the shift as a function of rotational longitude 
on 2005 FY9 for a particular CH 4 -ice ban d. Such a techn i que c ould provide a measure of 



inhomogeneity on the surface of 2005 FY9. iGrundy fc Finkl (119961 ) made such measurements 



for the 8897 A band on Pluto and found that the blueshift varied from to 10 A over the 
surface. 

On Pluto, it appears N 2 is sufficiently mobile to form textures that allow incoming 
photons to travel several centimeters through the very transparent, polycrystalline N 2 -ice 
before the photons are scattered off inclusions or grain boundaries. The long pathlengths 
make it possible for a small amount of CH 4 relative to N 2 to give deep CH 4 absorption bands. 
The mobility of CO is only a little less than N 2 , and the mobility of Ar is only a little less 
than CO. Therefore, it is possible that CO-ice or Ar-ice could provide the host matrix for 
long optical pathlengths. It is also possible the host matrix is a combination of N 2 , CO, and 
Ar. A suite of laboratory experiments that measure the shifts of CH 4 absorption maxima 
for different concentrations of CH 4 relative to N 2 , CO, and Ar would be highly valuable for 
the future interpretation of observational data. 
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Our spectrum of 2003 EL61 is remarkably featureless. There is a hint of a band at 
5773 A possibly due to 2 -ice. Future work will attempt to double the signal precision 
of the continuum to at least 300, which is the level that appears to be required for any 
possible detection of the weaker 02-ice band at 6275 A. The detection of 02-ice in significant 
quantities on KBOs could be of considerable interest in the distant future as a source of 
spacecraft fuel. 
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Fig. 1. — Our normalized reflectance spectrum of 2005 FY9 (bottom). Tick marks indicate 
the position o f previously reported CH 4 -ice bands, see Table 1. The middle spectrum is als o 
of 2005 FY9 jLicandro et alibood ). The top spectrum is of Pluto (torundv fc Finkl ll996l). 
All spectra are normalized to 1 at 6500 A. The iLicandro et al.l (120061 ) and Pluto spectra are 
offset by 0.4 and 0.8, respectively. Clearly, 2005 FY9 has deeper CH4-ice bands than Pluto. 
Previously unreported CH 4 -ice bands appear near 5400, 5800, and 6200 A in our spectrum 
of 2005 FY9. 
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Fig. 2. — Approximate albedo spectrum of 2005 FY9 (black line) and a Hapke model with 
a reddening agent and a CH 4 -ice grain size of 0.55 cm (red line). We assume an albedo at 
6500 A of 0.8. The model fits the core of the 8897 A band, but has too little absorption at 
all other bands. There are no laboratory absorption coefficients to fit CH 4 bands with A < 
7000 A. 
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Fig. 3. — Approximate albedo spectrum of 2005 FY9 (black line) and a Hapke model with 
a reddening agent and two CH4-ice grain sizes, 6 cm and 0.1 cm (red line). We assume an 
albedo at 6500 A of 0.8. Grains sizes larger than a cm probably provide a measure of the 
spacing between fractures or voids in the surface rather than indicating a surface covered 
with golf-ball sized particles. There are no laboratory absorption coefficients to fit CH 4 
bands with A < 7000 A 
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Fig. 4. — The minimum of the 7296 A CH 4-ice band in 2005 FY9 (black line) is blueshifted 
relative to the minimum of the pure CH4-ice Hapke model (red line). The other bands in 2005 
FY9 show a similar shift. The shift may be due to the presence of another ice component 
on the surface of 2005 FY9. 
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Fig. 5. — Result of cross-correlation experiment between the spectra of 2005 FY9 and the 
pure CH 4 -ice Hapke model in Figure 3. The difference between the two spectra reach a 
minimum at a 3.25 ± 2.25 A blueshift of the data relative to the model. 
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Fig. 6. — Normalized refle ctance spectrum of 2003 EL61 (black line) and the spectral ratio 
of Ganymede to Callisto (jCalvin fc Spencerlll997l ) showing 2 -ice bands at 5773 and 6275 
A (red line). There is a tantalizing dip in the spectrum of 2003 EL61 at the position of the 
5773 A band; however, a much higher continuum signal precision is necessary to test for the 
presence of the weaker 6275 A band. 
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Table 1. Laboratory Wavelengths and Frequencies of CH 4 -Ice Absorption Maxima 



Transition 


Wavelength (A) 


Wavenumber (cm 1 ) 


3^i + 4i/4 


7296 


13706 


3u 3 + 3u4 


7862 


12719 


3vi + 3vA 


7993 


12511 


4i/ 3 


8415 


11884 


vi + 3u 3 


8442 


11846 


3u 3 + 2U4. 


8691 


11506 


2ui + u 3 + 2u 4 


8897 


11240 


3u± + 2u4 


8968 


11151 


2v>, + 4U4. 


9019 


11088 



a Pure CH 4 Ice I at 30 K (Grundy et al. 2002) 



